Detection of Accelerator-Produced Neutrinos at a Distance of 250 km 
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The KEK to Kamioka long-baseline neutrino experiment (K2K) has begun its investigation of 
neutrino oscillations suggested by atmospheric neutrino observations. Twenty-eight neutrino events 
have been detected in coincidence with the expected arrival time of the beam in the 22.5 kt fiducial 
volume of Super-Kamiokande, the far detector at 250 km distance. The expectation is 37.8l3 g, 
derived using measurements of neutrino interactions in a near detector and extrapolation using a 
beam simulation validated by a measurement of pion kinematics after production and focusing. The 
background is of order 10""^ events. 

PACS numbers: 14.60.Pq, 13.15.-hg, 23.40.Bw, 95.55. Vj 
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The Standard Model of electroweak interactions has 
been tested with exceptional precision. The model, how- 
ever, does not address the question of the origin of gen- 
erations and their mixing. The existence of neutrino 
oscillations implies that neutrinos are massive and that 
lepton flavors are not conserved quantum numbers. It 
enables us to study relations between mass and flavor 
eigenstates in the lepton sector. Experiments on at- 
mospheric neutrinos have found a significant deficit in 
the flux of which have travelled an earth scale dis- 
tance [^||. The interpretation of these results provides 
not only strong evidence of Ut (or j/^ — > I'stcriic) 

oscillations |^,^ but also evidence for different mixing 
properties in the lepton and quark sectors. For a neu- 
trino energy Ey{GeY) and a distance from the source 
L(km), the oscillation probability can be written in 
terms of the mixing angle 9 and the difference of the 
squared masses Am^ (eV^) in two flavor approximation 
as P{v^ v^) = sin^ 2B aiw^ [1.27 /^m^L/E^). 

The KEK to Kamioka long-baseline neutrino exper- 
iment (K2K) is the first accelerator-based experiment 
with hundreds of km neutrino path length. The intense, 
nearly pure neutrino beam (98.2% v^, 1.3% v^, and 0.5% 
Vfj) has an average L/E^ « 200(L — 250km, {E^) « 
1.3GeV). The neutrino beam properties are measured 
just after production, and the kinematics of parent pions 
are measured in situ to extrapolate the measurements 
at the near detector to the expectation at the far detec- 
tor. K2K focuses on the study of the existence of neu- 
trino oscillations in disappearance that is observed in 
atmospheric neutrinos, and on the search for — > v,, 
oscillation with well understood flux and neutrino com- 
position in the Am^ > 2 x 10^'^eV^ region. In this let- 
ter, the first results on the event rates with 100 days of 
data taken from June 1999 to June 2000 corresponding 
to 2.29 X 10^^ protons on target, and the performance 
of the critical components of the long-baseline neutrino 
experiment are described. 

The K2K neutrino beam |6[0] is a horn-focused wide 
band beam. The primary beam for K2K is 12 GeV ki- 
netic energy protons from the KEK proton-synchrotron 
(KEK-PS) [§. Every 2.2 s, approximately 6 x 10^^ pj-^. 
tons in nine bunches are fast-extracted in a single turn, 
making a 1.1 /iS beam spill. These protons are focused 
onto a 30 mm diameter, 66 cm long aluminum target 
which is a current carrying element in the first of a pair 
of horn magnets operating at 250 kA This design 
maximizes the efficiency of the toroidal magnetic field to 
focus positive pions produced in proton-aluminum inter- 
actions, while sweeping out negative secondary particles. 

Downstream of the horn system, before the 200 m long 
decay volume where the pions decay to and muons, 
a gas-Cherenkov detector (PIMON) |l^] is occasionally 
put in the beam to measure the kinematics of the pions 
after their production and subsequent focusing. After the 
decay volume, an iron and concrete beam dump stops 



essentially all charged particles except muons with an 
energy greater than 5.5 GeV. Downstream of the dump 
there is a "muon monitor" (MUMON) ||l^ consisting of a 
segmented ionization chamber and an array of silicon pad 
detectors. This monitors the residual muons spill-by-spill 
to check beam centering and muon yield. The ionization 
chamber consists of 5 cm wide strips covering roughly 2x 
2m^ area with separate planes for horizontal and vertical 
read-out. The 26 silicon pads, each of ~ lOcm^ transverse 
area, are distributed through a 3 x 3m^ area. 

Neutrino interactions near the production site are mea- 
sured by a set of detectors with complementary abilities. 
The detectors are located approximately 300 m from the 
pion production target with approximately 70 m of earth 
eliminating virtually all prompt beam products other 
than neutrinos. A one kiloton water Cherenkov detector 
(Ikt) uses the same technology and analysis algorithms 
as the far detector, Super-Kamiokande (SK). It has 680 
20" photomultiplier tubes (PMTs) on a 70 cm grid lining 
a 8.6m diameter, 8.6m high cylinder. The PMTs them- 
selves and their arrangement are the same as in SK, giv- 
ing the same fractional coverage by photo-cathode (40%). 
A scintillating fiber detector (SciFi) with a 6 ton 
water target has tracking capability, and allows discrim- 
ination between different types of interactions such as 
quasi-elastic or inelastic. Downstream of the SciFi there 
is a lead glass array for tagging electromagnetic showers. 
A muon range detector (MRD) |jl^, measures the energy, 
angle, and production point of muons from charged cur- 
rent (CC) interactions. It covers 7.6x7.6m^ transverse 
area with four 10-cm iron plates, followed by eight 20-cm 
iron plates, all interleaved with drift tubes. The total 
mass is 915 tons. The 6632 drift tubes each has 5x7cm^ 
cross section, and are arranged in horizontal and vertical 
read-out planes. 

The far detector for K2K is the SK detector located 
in the Kamioka Observatory, Institute for Cosmic Ray 
Research (ICRR), University of Tokyo, which has been 
taking data since 1996. Its performance and results have 
been documented in the literature Q . Event selection for 
this detector uses timing synchronization with the KEK- 
PS via the Global Positioning System (GPS) 

The beam- line was aligned by GPS position survey 
The precision of this survey for the line from the target 
to the far detector is better than 0.01 mr and the con- 
struction precision for the near site alignment is better 
than 0.1 mr. The predicted neutrino spectrum at 250 
km is approximately the same over nearly 1 km, giving 
a required accuracy of 3mr for the pointing of the beam. 

The steering and the monitoring of the beam are 
carried out based on the MUMON measurement. The 
and muons in the decay volume originate from the 
same pion decays, so the measured center position of the 
muon profile is correlated with the i/^ beam direction. 
The r.m.s. muon profile width is ^ 90cm for both the 
horizontal and vertical directions. These measurements 
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are used in tuning the beam direction at the start of ev- 
ery beam period, usually once per month. Since there is 
a magnetic field in and around the target, the direction 
of the secondary pions, which determines the direction 
of decay and muons, is sensitive to the position of 
the primary proton beam at the target. A 1 mm dis- 
placement of protons at the target face causes roughly 
0.5 mr deflection to the opposite direction for muons ob- 
served by the MUMON Q. During beam tuning, the 
primary proton position at the target is adjusted so that 
the muon profile is centered on the SK direction within 
0.1 mr. The center of the muon profile is measured to 
be stable spill-by-spill within ±lmr in each dimension 
throughout the whole experimental period. The stability 
of the muon yield is directly related to the stability of 
the horn-focused beam. The yield normalized to pro- 
ton intensity measured by a current transformer is stable 
spill-by-spill within a measurement uncertainty of 2%. 

The characteristics and stability of the beam itself 
are directly monitored at the near site by the MRD us- 
ing interactions with iron. The large transverse area 
of MRD makes it possible to measure the beam direction 
and width (profile), and the large mass makes it possi- 
ble to measure the time stability of the event rate, 
profile and spectrum. The location of the center of the 
profile gives the beam direction. For the profile measure- 
ment, the starting point of a reconstructed muon track 
is regarded as the vertex of the f^-iron interaction. The 
data reduction for the profile measurement requires: a) 
Only tracks around the time of the beam spill are ac- 
cepted, b) Tracks with a common vertex are taken as 
a single event with the longest track assumed to be a 
muon. c) Muons entering or exiting the detector are re- 
jected, d) Any muons with reconstructed energy lower 
than 0.5 GeV or higher than 2.5 GeV are rejected in 
order to avoid regions of small acceptance. The vertex 
distribution is corrected for geometrical acceptance to 
obtain the beam profile. As shown in Fig. |^(a), the hor- 
izontal profile is well centered on the SK direction. The 
profile width is well reproduced by the beam simulation 
which is described below. The profile center is plotted 
as a function of time in Fig. |l](b). The i/^ direction has 
been stable within ±1 mr throughout all data-taking pe- 
riods. The vertical profile is also well centered on the SK 
direction and stable. For the event rate and kinematic 
distribution analyses no energy cut is made, but only 
events with vertices inside the fiducial volume defined by 
a cylinder of 3m radius in the upstream 9 iron plates are 
accepted. The measured rate of j/^-iron events is typi- 
cally 0.05 events/spill and can be monitored on a daily 
basis with good statistics. The event rate normalized to 
proton intensity is stable within the statistical error of 
the MRD measurement. The muon energy and angu- 
lar distributions are also continuously monitored. They 
show no change as a function of time, implying the v^^ 
energy spectrum is stable throughout this period. 
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FIG. 1. (a) Horizontal profile of the z/^ beam measured by 
MRD in a one month period (points with error bars), and 
simulated (boxes; size corresponds to error). Errors include 
statistics and uncertainty of the acceptance correction. The 
normalization is by the area under the histograms, (b) Sta- 
bility of the Vf_i beam direction. The fitted center of the hor- 
izontal profile is plotted for five day periods, with error bars 
from the fits. The dashed lines show rtlmr directions to SK. 



Comparison between expected and observed numbers 
of events at the far site is done with this knowledge of 
measured beam stability. To predict the beam char- 
acteristics at the far site, a normalization measurement 
at the near site and the extrapolation of the information 
from near to far are necessary. For the rate normaliza- 
tion, the Ikt is used so that any detector or analysis bias 
is suppressed. For the extrapolation, the beam simula- 
tion is used, which is validated by the PIMON measure- 
ment as described below. This simulation is based on 



GEANT 1 15 1 with detailed description of materials and 
magnetic fields in the target region and decay volume. It 
uses as input a measurement of the primary beam profile 
at the target. Primary proton interactions on aluminum 
are modeled with a parameterization of hadron produc- 
tion data Other hadronic interactions are treated 
by GEANT-CALOR 0. 

Once the kinematic distribution of pions after produc- 
tion and focusing is known, it is possible to predict the 
spectrum at any distance from the source. The pion 
momentum and angular distribution is measured by the 
PIMON. Cherenkov photons generated by pions in the 
detector are focused by a pie-shaped spherical mirror to 
an array of PMTs of 8 mm effective diameter in its focal 
plane. Photons emitted by particles with the same veloc- 
ity and incident angle with respect to the mirror arrive at 
the same position in the focal plane. This is independent 
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of photon incident position with respect to the mirror. 
The photon distribution on the PMT array is a superpo- 
sition of shoes of the Cherenkov rings from particles of 
various velocities and angles. Measurements are made at 
seven indices of refraction, controlled by gas pressure, to 
give additional information on the velocity distribution 
of the particles. The indices of refraction are chosen so 
that the corresponding pion momentum interval is 400 
MeV/c. The pion two-dimensional distribution of mo- 
mentum versus angle is derived by unfolding the photon 
distribution data with various indices of refraction. The 
binning of extracted pion momentum and angular distri- 
butions are 1 GeV/c and 10 mr, respectively. To avoid 
background from 12 GeV primary protons, the index of 
refraction is adjusted below the Cherenkov threshold of 
12 GeV protons (1 — /? w2.6 x 10"'^). As a consequence, 
analysis is done for > 2GeV/c, giving the energy 
spectral shape above 1 GeV. Figure ^ shows the inferred 
i/^ energy spectral shape at the near site and the far to 
near flux ratio along with the beam simulation result. 
The beam simulation is well validated by the PIMON 
measurement without any tuning. 
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FIG. 2. The top figure shows the i/^ energy spectral shape 
at the near site and the bottom figure shows the far to near 
1/^ flux ratio. The histograms are from the beam simulation 
results. The data points are derived from the PIMON mea- 
surement. For the top figure, normalization is done by area 
above 1 GeV. The vertical error bars for the data points refiect 
the total uncertainty. 

The Vfj. interaction rate at the near site is measured in 
the Ikt by detecting Cherenkov light emitted from pro- 
duced charged particles. The vertices and directions of 
Cherenkov rings are reconstructed with the same meth- 
ods as in SK [^. In the Ikt, however, an analog sum 
of signals from all PMTs is recorded by an FADC to se- 
lect beam spills with only one event. The definition of 
events per spill is the number of peaks of the FADC signal 
with > 1000 collected photo-electrons (p.e.) (« lOOMeV 



deposited energy). The average number of events ob- 
served in the full detector is about 0.2/spill including 
entering background events due to cosmic rays or in- 
duced muons from upstream (about 1.5% of events in the 
fiducial volume defined below). Thus >2 events occur in 
about 10% of those spills with > 1 event. Neutrino inter- 
actions are selected by requiring: a) There is no detector 
activity in the 1.2/xs preceding the beam spill, b) Only 
a single event is observed in the FADC peak search for 
that spill, c) The reconstructed vertex is inside the 25 t 
fiducial volume defined by a 2 m radius, 2 m long cylinder 
along the beam axis, in the upstream side of the detector. 
The detection efficiency of the Ikt detector for detecting 
CC interactions is 87% and for neutral current inelastic 
(NCinei) interactions it is 55%. The expected ratio of CC 
to NCinci gives an overall efficiency of 72%. The main 
source of inefficiency is the 1000 p.e. threshold of the 
peak search in the FADC signal. The neutrino interac- 
tion simulator used for efficiency calculations is the same 
as that used in all SK analyses. Finally the number of 
events in Ikt is corrected for spills with multiple events. 
The average j/^ event rate per proton hitting the target 
is 3.2 X 10~^^. Correcting for efficiencies, relative target 
masses, and detector live times, the expected signal 
at the far site is estimated by applying the extrapolation 
from the experimentally validated beam simulation. 

Accelerator-produced neutrino interactions at the far 
site, SK, are selected by comparing two Universal Time 
Coordinated (UTC) time stamps from the GPS system, 
Tkek for the KEK-PS beam spill start time and TgK for 
the SK trigger time. The time difference between two 
UTC time stamps, A(T) = Tsk - Jkek - TOF where 
TOF is the time of flight of a neutrino, should be dis- 
tributed around the interval from to 1.1/is to match the 
width of the beam spill of the KEK-PS. Since the mea- 
sured uncertainty of the synchronization accuracy for the 
two sites is < 200ns, beam-induced events are selected in 
a 1.5^s window. Data reduction similar to that used in 
atmospheric neutrino analyses at SK is applied to 
select fully contained (FC) neutrino interactions. The 
criteria are: a) There is no detector activity within 30 
fis before the event, b) The total collected p.e. in a 
300 ns time window is > 200 (w 20MeV deposited en- 
ergy), c) The number of PMTs in the largest hit cluster 
in the outer-detector is < 10. d) The deposited energy is 
> 30 MeV. Finally, a fiducial cut is applied accepting only 
events with fitted vertices inside the same 22.5kt volume 
used for SK atmospheric neutrino analysis. The detec- 
tion efficiency of SK is 93% for CC interactions and 68% 
for NCinci interactions, for a total of 79%. Similarly to 
the Ikt, the inefficiency is mainly due to the energy cut. 
Figure |^ shows the A(r) distribution at various stages 
of the reduction. A clear peak in time with the neutrino 
beam from the KEK-PS is observed in the analysis time 
window. Twenty-eight FC events are observed in the 
fiducial volume. The arrival rate of neutrinos observed 
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by the far detector is consistent with constant within sta- 
tistical fluctuation. The 1.5/xs selection window gives an 
expected background from atmospheric neutrino interac- 
tions of order 10^^ events. 



TABLE I. Summary of the observed and expected numbers 
of FC events in the fiducial volume at the far site. 
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3. Far site timing analysis; (top) A(T) distribution 
over a zt500/is window for events remaining after cuts a) (un- 
filled histogram) and b) (hatched), and after all cuts except 
timing (shaded). The cuts are described in the text, (bot- 
tom) A(T) for a ±5^s window for FC events with vertices in 
the fiducial volume. 

The systematic uncertainty of the Ikt measurement is 
5%, for which the leading terms are due to vertex fitting 
and its effect on the fiducial cut (4%), and the treatment 
of multiple events in a spill (3%). Other sources such as 
energy scale uncertainty are relatively small. The statis- 
tical uncertainty of the Ikt measurement is < 1%. The 
systematic uncertainty of the extrapolation from the near 
site measurement is estimated to be iy% based on the 
PIMON measurement uncertainty and beam simulation 
uncertainty for low energy neutrinos. The systematic 
uncertainty in the SK measurement is 3%, mainly due 
to the fiducial cut. The systematic uncertainty term 
coming from uncertainty in the neutrino energy spec- 
trum and cross section is small due to cancellations. The 
quadrature sum of all known uncertainties for the far 
site event rate prediction is tio%- The resulting expec- 



tation is 37.8_3 g events in the absence of neutrino os- 
cillations. The observed and expected numbers for var- 
ious categories are summarized in Table |. The number 
of Cherenkov rings and particle identification are recon- 
structed by the same algorithms as those used at SK , 
and the event category definitions are also the same. The 
expected number of events for analyses of CC interactions 
in the other near detectors are 41.0l|^ for MRD ( non 
target) and 37. 2'^^^'^ for SciFi (water target) ||l^], each of 
which is consistent with the value based on Ikt analysis. 

The crucial principles of K2K have been proven to 
work: The beam direction has been continuously mon- 
itored and controlled within ±lmr for several months 
duration by controlling the proton position on target 



Event Category 
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Single Ring /i-like 


14 


20.9 


Single Ring e-like 
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Multi Ring 


13 


14.9 


TOTAL 


28 


37.8^3:^8 



with millimeter accuracy and operating the horn mag- 
net stably. Pion kinematics have successfully been mea- 
sured in situ allowing prediction of the neutrino beam 
at the distance of 250 km with six to seven percent ac- 
curacy. Events at the far site have been identified by 
means of GPS timing, reducing backgrounds to a negli- 
gible level. Twenty-eight FC neutrino events have been 
observed where 37.8I3 g are expected. The experiment 
expects to accumulate 10^" protons on target, provid- 
ing sufficient statistics to study neutrino oscillations by 
spectral analysis for disappearance. 
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